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CJ Abstract 

The fluctuations of two-dimensional extended objects {membranes) 
^1 is a rich and exciting field with many solid results and a wide range of 

OO ' open issues. We review the distinct universality classes of membranes, 

■ determined by the local order, and the associated phase diagrams. Af- 
i ter a discussion of several physical examples of membranes we turn to 
I the physics of crystalline (or polymerized) membranes in which the 

individual monomers are rigidly bound. We discuss the phase dia- 
gram with particular attention to the dependence on the degree of 
I self-avoidance and anisotropy. In each case we review and discuss an- 

alytic, numerical and experimental predictions of critical exponents 
. and other key observables. Particular emphasis is given to the results 

', obtained from the renormalization group e-expansion. The resulting 

renormalization group flows and fixed points are illustrated graphi- 
cally. The full technical details necessary to perform actual calcula- 
tions are presented in the Appendices. We then turn to a discussion of 
^ ' the role of topological defects whose liberation leads to the hexatic and 

■ fluid universality classes. We finish with conclusions and a discussion 
of promising open directions for the future. 
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1 Introduction 



The statistical mechanics of one- dimensional structures (polymers) is fasci- 
nating and has proved to be fruitful from the fundamental and applied points 
of view [jl], . The key reasons for this success lie in the notion of universality 
and the relative simplicity of one-dimensional geometry. Many features of 
the long- wavelength behavior of polymers are independent of the detailed 
physical and chemical nature of the monomers that constitute the polymer 
building blocks and their bonding into macromolecules. These microscopic 
details simply wash out in the thermodynamic limit of large systems and 
allow predictions of critical exponents that should apply to a wide class of 
microscopically distinct polymeric systems. Polymers are also sufficiently 
simple that considerable analytic and numerical progress has been possible. 
Their statistical mechanics is essentially that of ensembles of various classes 
of random walks in some (i- dimensional bulk or embedding space. 

A natural extension of these systems is to intrinsic two-dimensional struc- 
tures which we may call generically call membranes. The statistical me- 
chanics of these random surfaces is far more complex than that of polymers 
because two-dimensional geometry is far richer than the very restricted ge- 
ometry of lines. Even planar two-dimensional — monolayers — are complex, 
as evidenced by the KTNHY [|], H, ^ theory of defect-mediated melting 
of monolayers with two distinct continuous phase transitions separating an 
intermediate hexatic phase, characterized by quasi-long-range bond orien- 
tational order, from both a low-temperature crystalline phase and a high- 
temperature fluid phase. But full-fledged membranes are subject also to 
shape fluctuations and their macroscopic behavior is determined by a sub- 
tle interplay between their particular microscopic order and the entropy of 
shape and elastic deformations. For membranes, unlike polymers, distinct 
types of microscopic order (crystalline, hexatic, fluid) will lead to distinct 
long-wavelength behavior and consequently a rich set of universality classes. 

Flexible membranes are an important member of the enormous class of 
soft condensed matter systems 0, ||, |^ , those which respond easily to exter- 
nal forces. Their physical properties are to a considerable extent dominated 
by the entropy of thermal fluctuations. 

In this review we will describe some of the presently understood behavior 
of crystalline (fixed-connectivity), hexatic and fiuid membranes, including 
the relevance of self-avoidance, intrinsic anisotropy and topological defects. 
Emphasis will be given to the role of the renormalization group in elucidating 
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the critical behavior of membranes. The polymer pastures may be lovely but 
a dazzling world awaits those who wander into the membrane meadows. 

The outline of the review is the following. In sec. ^ we describe a vari- 
ety of important physical examples of membranes, with representatives from 
the key universality classes. In sec. ^ we introduce basic notions from the 
renormalization group and some formalism that we will use in the rest of the 
review. In sec. H we review the phase structure of crystalline membranes for 
both phantom and self-avoiding membranes, including a thorough discussion 
of the fixed-point structure, RG flows and critical exponents of each global 
phase. In sec. ^we turn to the same issues for intrinsically anisotropic mem- 
branes, with the new feature of the tubular phase. In sec. ^ we address the 
consequences of allowing for membrane defects, leading to a discussion of the 
hexatic membrane universality class. We end with a brief discussion of fluid 
membranes in sec. |^ and conclusions. 



2 Physical examples of membranes 

There are many concrete realizations of membranes in nature, which greatly 
enhances the significance of their study. Crystalline membranes, sometimes 
termed tethered or polymerized membranes, are the natural generalization 
of linear polymer chains to intrinsically two-dimensional structures. They 
possess in-plane elastic moduli as well as bending rigidity and are character- 
ized by broken translational invariance in the plane and fixed connectivity 
resulting from relatively strong bonding. Geometrically speaking they have 
a preferred two-dimensional metric. Let's look at some of the examples. 
One can polymerize suitable chiral oligomeric precursors to form molecular 
sheets This approach is based directly on the idea of creating an intrinsi- 
cally two-dimensional polymer. Alternatively one can permanently cross-link 
fluid-like Langmuir-Blodgett films or amphiphilic bilayers by adding certain 
functional groups to the hydrocarbon tails and/or the polar heads |jTT|, as 
shown schematically in Fig. |I]. 

The cytoskeletons of cell membranes are beautiful and naturally occur- 
ring crystalline membranes that are essential to cell membrane stability and 
functionahty. The simplest and most thoroughly studied example is the cy- 
toskeleton of mammalian erythrocytes (red blood cells). The human body 
has roughly 5 x 10^'^ red blood cells. The red blood cell cytoskeleton is a 
fishnet-like network of triangular plaquettes formed primarily by the proteins 
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Polymerized membrane 



Figure 1: The polymerization of fluid- like membrane to a crystalline mem- 
brane. 



spectrin and actin. The links of the mesh are spectrin tetramers (of length 
approximately 200 nm) and the nodes are short actin filaments (of length 
37 nm and typically 13 actin monomers long) |13|, |T^, as seen in Fig.0 and 
Fig.^. There are roughly 70,000 triangular plaquettes in the mesh altogether 
and the cytoskeleton as a whole is bound by ankyrin and other proteins to 
the cytoplasmic side of the fluid phospholipid bilayer which constitutes the 
other key component of the red blood cell membrane. 

There are also inorganic realizations of crystalline membranes. Graphitic 
oxide (GO) membranes are micron size sheets of solid carbon with thicknesses 
on the order of lOA, formed by exfoliating carbon with a strong oxidizing 
agent. Their structure in an aqueous suspension has been examined by sev- 
eral groups []r^ |TB|, ^]. Metal dichalcogenides such as M0S2 have also been 
observed to form rag- like sheets [^. Finally similar structures occur in the 
large sheet molecules, shown in Fig.^, believed to be an ingredient in glassy 
B2O3. 

In contrast to crystalline membranes, fluid membranes are characterized 
by vanishing shear modulus and dynamical connectivity. They exhibit signif- 
icant shape fluctuations controlled by an effective bending rigidity parameter. 

A rich source of physical realizations of fluid membranes is found in am- 
phiphilic systems pT] , |22| , [23|] . Amphiphiles are molecules with a two-fold 
character - one part is hydrophobic and another part hydrophilic. The clas- 
sic examples are lipid molecules, such as phospholipids, which have polar or 
ionic head groups (the hydrophilic component) and hydrocarbon tails (the 
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Figure 2: An electron micrograph of a region of the erythrocyte cytoskele- 
ton. The skeleton is negatively stained (magnification 365,000) and has been 
artificially spread to a surface area nine to ten times as great as in the native 
membrane |]T5| . 





Figure 3: An extended view of the crystalline spectrin/actin network which 
forms the cytoskeleton of the red blood cell membrane |jl6| . 



Figure 4: The sheet molecule -B2O3 
WATER 




Figure 5: Schematic of experimental procedure to make a black membrane. 

hydrophobic component). Such systems are observed to self-assemble into 
a bewildering array of ordered structures, such as monolayers, planar (see 
Fig.^ and spherical bilayers (vesicles or liposomes) (see Fig. ||) as well as 
lamellar, hexagonal and bicontinuous phases [Q. In each case the basic 
ingredients are thin and highly flexible surfaces of amphiphiles. The lipid 
bilayer of cell membranes may itself be viewed as a fluid membrane with 
considerable disorder in the form of membrane proteins (both peripheral and 
integral) and with, generally, an attached crystalline cytoskeleton, such as 
the spectrin/actin mesh discussed above. 

A complete understanding of these biological membranes will require a 
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WATER 



WATER 



25 nm 

Figure 6: The structure of a liposome with its pure hpid spherical bilayer 

thorough understanding of each of its components (fluid and crystalline) 
followed by the challenging problem of the coupled system with thermal fluc- 
tuations, self-avoidance, potential anisotropy and disorder. The full system 
is currently beyond the scope of analytic and numerical methods but there 
has been considerable progress in the last fifteen years. 

Related examples of fluid membranes arise when the surface tension be- 
tween two normally immiscible substances, such as oil and water, is sig- 
nificantly lowered by the surface action of amphiphiles (surfactants), which 
preferentially orient with their polar heads in water and their hydrocarbon 
tails in oil. For some range of amphiphile concentration both phases can 
span the system, leading to a bicontinuous complex fluid known as a mi- 
croemulsion. The oil-water interface of a microemulsion is a rather unruly 
fluid surface with strong thermal fluctuations |^ (see Fig.|^). 

The structures formed by membrane/polymer complexes are of consider- 
able current theoretical, experimental and medical interest. To be specific 
it has recently been found that mixtures of cationic liposomes (positively 
charged vesicles) and linear DNA chains spontaneously self-assemble into a 
coupled two-dimensional smectic phase of DNA chains embedded between 



lamellar lipid bilayers |26, 27|. For the appropriate regime of lipid concen- 
tration the same system can also form an inverted hexagonal phase with 
the DNA encapsulated by cylindrical columns of liposomes[^] (see Fig.^. 



In both these structures the liposomes may act as non-viral carriers (vec- 
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Figure 7: The structure of a microemulsion formed by the addition of sur- 
factant to an oil-water mixture. 

METAL 




Figure 8: Metal-coated fluid microcylinders (tubules) formed by chiral lipids. 
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tors) for DNA with many potentially important applications in gene ther- 
apy p9|. Liposomes themselves have long been studied and utilized in the 
pharmaceutical industry as drug carriers BO]. On the materials science side 




the self-assembling ability of membranes is being exploited to fabricate mi- 
crostructures for advanced material development. One beautiful example is 
the use of chiral-lipid based fluid microcylinders (tubules) as a template for 
metallization. The resultant hollow metal needles may be half a micron in di- 
ameter and as much as a millimeter in length ||31], ^ , as illustrated in Fig.^. 
They have potential applications as, for example, cathodes for vacuum field 



emission and microvials for controlled release 31 



3 The Renormalization Group 

The Renormalization Group (RG) has provided an extremely general frame- 
work that has unified whole areas of physics and chemistry . It is beyond 
the scope of this review to discuss the RG formalism in detail but there is an 
ample literature to which we refer the reader (see the articles in this issue). 
It is the goal of this review to apply the RG framework to the statistical me- 
chanics of membranes, and for this reason we briefly emphasize and review 
some well known aspects of the RG and its related ^-expansion. 
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The RG formalism elegantly shows that the large distance properties 
(or equivalently low p-limit) of different models are actually governed by 
the properties of the corresponding Fixed Point (FP). In this way one can 
compute observables in a variety of models, such as a molecular dynamics 
simulation or a continuum Landau phenomenological approach, and obtain 
the same long wavelength result. The main idea is to encode the effects of 
the short-distance degrees of freedom in redefined couplings. A practical way 
to implement such a program is the Renormalization Group Transformation 
(RGT), which provides an explicit prescription for integrating out all the 
high p-modes of the theory. One obtains the large-distance universal term of 
any model by applying a very large (oo to be rigorous) number of RGTs. 

The previous approach is very general and simple but presents the tech- 
nical problem of the proliferation in the number of operators generated along 
the RG flow. There are established techniques to control this expansion, one 
of the most successful ones being the £-expansion. The £-expansion may 
also performed via a field theoretical approach using Feynman diagrams and 
dimensional regularization within a minimal subtraction scheme, which we 
briefly discuss below. Whereas it is true that this technique is rather ab- 
stract and intuitively not very close to the physics of the model, we find it 
computationally much simpler. 

Generally we describe a particular model by several fields {(f), x, • • •} and 
we construct the Landau free energy by including all terms compatible with 
the symmetries and introducing new couplings [u, v, ■ ■ ■) for each term. The 
Landau free energy may be considered in arbitrary dimension d, and then, one 
usually finds a Gaussian FP (quadratic in the fields) which is infrared stable 
above a critical dimension (du)- Below du there are one or several couplings 
that define relevant directions. One then computes all physical quantities as 
a function of e = djj — d, that is, as perturbations of the Gaussian theory. 

In the field theory approach, we introduce a renormalization constant for 
each field {Z^, Z^, ■ ■ ■) and a renormalization constant {Zu, Zy, ■ ■ •) for each 
relevant direction below du- If the model has symmetries, there are some 
relations among observables (Ward identities) and some of these renormal- 
ization constants may be related. This not only reduces their number but 
also has the added bonus of providing cross-checks in practical calculations. 
Within dimensional regularization, the infinities of the Feynman diagrams 
appear as poles in £, which encode the short-distance details of the model. If 
we use these new constants (^'s) to absorb the poles in e, thereby producing 
a complete set of finite Green's functions, we have succeeded in carrying out 
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the RG program of including the appropriate short-distance information in 
redefined couphngs and fields. This particular prescription of absorbing only 
the poles in e in the Z's is called the Minimal Subtraction Scheme (MS), and 
it considerably simplifies practical calculations. 

As a concrete example, we consider the theory of a single scalar field 
(f) with two independent coupling constants. The one-particle irreducible 
Green's function has the form 



(k,; UR, VR, M) = Z^/'r^(k,; u, v; ^ 



1, 



(1) 



where the function on the left depends on a new parameter M, which is 
unavoidably introduced in eliminating the poles in e. The associated correla- 
tor also depends on redefined couplings ur and vr. The rhs depends on the 
poles in e, but its only dependence on M arises through Z^. This observation 
allows one to write 



M 



R 



_d_ /^-iV/2p(iV) 

dM V ^ 



(2) 



dM 



du 



R 



_d_ 

dvf 



where 



0, 



UR = M-'F{Z^,Z^,---\Zu)u 



(3u{ur,Vr) 



M 



dM 



VR = M-'F{Z^,Z^r--\Zv)v 
dM, 



70 




The /5-functions control the running of the coupling by 



M 



du 



R 



dv 



f3u{uR,VR) , = f3^{uR,VR) 



dM 



dM 



(3) 



(4) 



The existence of a FP, at which couplings cease to flow, requires (3{u\, f |j) = 
for all /5-functions of the model. Those are the most important aspects of 
the RG we wanted to review. In Appendix |B] we derive more appropriate 
expressions of the RG-functions for practical convenience. For a detailed 
exposition of the e-expansion within the field theory framework we refer to 



the excellent book by Amit[34| 
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4 Crystalline Membranes 



A crystalline membrane is a two dimensional fish-net structure with bonds 
(links) that never break - the connectivity of the monomers (nodes) is fixed. 
It is useful to keep the discussion general and consider D-dimensional objects 
embedded in c?-dimensional space. These are described by a d-dimensional 
vector r(x), with x the D-dimensional internal coordinates, as illustrated 
in Fig.p!0|. The case [d = 3,D = 2) corresponds to the physical crystalline 
membrane. 

Zt 




Figure 10: Representation of a membrane. 



To construct the Landau free energy of the model, one must recall that 
the free energy must be invariant under global translations, so the order 
parameter is given by derivatives of the embedding r, that is ta = with 
a = 1,---,D. This latter condition, together with the invariance under 



rotations (both in internal and bulk space), give a Landau free energy ||35 



Fir) 



^/rf^xciVV^-ny)) 



(5) 
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where higher order terms may be shown to be irrelevant at long wavelength, 
as discussed later. The physics in Eq.(^) depends on five parameters, 



• ft, bending rigidity : This is the coupling to the extrinsic curvature 
(the square of the Gaussian mean curvature). Since reparametrization 
invariance is broken for crystalline membranes, this term may be re- 
placed by its long- wavelength limit. For large and positive bending 
rigidities fiatter surfaces are favored. 

• t. u. V. elastic constants : These coefficients encode the microscopic elas- 
tic properties of the membrane. In a fiat phase, they may be related 
to the Lame coefficients of Landau elastic theory (see sect. |4.1.3| ). 

• b, Excluded volume or self-avoiding coupling : This is the coupling 
that imposes an energy penalty for the membrane to self- intersect. The 
case 6 = 0, i. e. no self- avoidance, corresponds to a phantom model. 

We generally expand r'(x) as 

f(x) = (Cx + u(x),/i(x)) , (6) 

with u the D-dimensional phonon in-plane modes, and h the d — D out-of- 
plane fiuctuations. If C = the model is in a rotationally invariant crumpled 
phase, where the typical surfaces have fractal dimension, and there is no 
real distinction between the in-plane phonons and out-of plane modes. For 



a pictorial view, see cases a) and b) in Fig .p. 

If C 7^ the membrane is fiat up to small fiuctuations and the full ro- 
tational symmetry is spontaneously broken. The fields h are the analog of 
the Goldstone bosons and they have different naive scaling properties than 
u. See Fig.^ for a visualization of a typical configuration in the fiat phase. 

We will begin by studying the phantom case first. This simplified model 
may even be relevant to physical systems since one can envision membranes 
that self-intersect (at least over some time scale). One can also view the 
model as a fascinating toy model for understanding the more physical self- 
avoiding case to be discussed later. Combined analytical and numerical stud- 
ies have yielded a thorough understanding of the phase diagram of phantom 
crystalline membranes. 
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Figure 11: Examples of a) crumpled phase, b) crumpling transition (crum- 
pled phase) and c) a flat phase. Results correspond to a numerical simulation 
of the phantom case |^ and gives a very intuitive physical picture of the 
different phases. 



4.1 Phantom 

The Phantom case corresponds to setting 6 = in the free energy Eq.(^): 



F(f) 



(7) 



The mean field effective potential, using the decomposition of Eq.(|]), 
becomes 

V{C) = DC{l: + {u + vD)e) , (8) 
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with solutions 

t 4{u+vD) ■ ^ U . 

There is, consequently, a fiat phase for t < and a crumpled phase for t > 0, 
separated by a crumpling transition at t = (see Fig.p!2D. 





Figure 12: Mean field solution for crystalline membranes. 



The actual phase diagram agrees qualitatively with the phase diagram of 
the model shown schematically in Fig.|13|. The crumpled phase is described 
by a line of equivalent FPs(GFP). There is a general hyper- surface, whose 
projection onto the n—t plane corresponds to a one-dimensional curve (CTH), 
which corresponds to the crumpling transition. Within the CTH there is an 
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infrared stable FP (CTFP) which describes the large distance properties of 
the crumpling transition. Finally, for large enough values of k and negative 
values of t, the system is in a flat phase described by the corresponding infra- 
red stable FP (FLFP) []. Although the precise phase diagram turns out to 



be slightly more complicated than the one depicted in Fig.[l^, the additional 
subtleties do not modify the general picture. 



CRUMPLED PHASE 



CTH 

(Crumping Transition) 
FLFP 




Figure 13: Schematic plot of the phase diagram for phantom membranes. 
GFP are the equivalent FPs describing the crumpled phase. The crumpling 
transition is described by the Crumpling transition critical line (CTH), which 
contains the Crumpling Transition FP (CTFP). The Flat phase is described 
by the (FLFP). 



The evidence for the phase diagram depicted in Fig.|l^ comes from com- 
bining the results of a variety of analytical and numerical calculations. We 
present in detail the results obtained from the e-expansion since they have 
wide applicability and allow a systematic calculation of the /5-function and 
the critical exponents. We also describe briefly results obtained with other 
approaches. 



The FLFP is actually a line of equivalent fixed points. 
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4.1.1 The crumpled phase 

In the crumpled phase, the free energy Eq.(^ for D > 2 simplifies to 

F{r) = ^ y d^x{dar)'^ + Irrelevant terms , (10) 

since the model is completely equivalent to a linear sigma model in D < 2 
dimensions having 0{d) symmetry, and therefore all derivative operators in r* 
are irrelevant by power counting. The parameter t labels equivalent Gaussian 
FPs, as depicted in Fig.|l^. In RG language, it defines a completely marginal 
direction. This is true provided the condition t > is satisfied. The large 
distance properties of this phase are described by simple Gaussian FPs and 
therefore the connected Green's function may be calculated exactly with 
result 

("^^^ ~ { loiTxl D t 2 (") 

The associated critical exponents may also be computed exactly. The Haus- 
dorff dimension dn-, or equivalently the size exponent v = D/dn, is given 
(for the membrane case D = 2) hj 

dn = oo {u = 0) ^ Rlr^logL . (12) 

The square of the radius of gyration Rq scales logarithmically with the mem- 
brane size L. This result is in complete agreement with numerical simulations 
of tethered membranes in the crumpled phase where the logarithmic behav- 
ior of the radius of gyration is accurately checked PPI, EO, ET], E3, HI, 



46t 0, 0. Reviews may be found in ||4J 



4.1.2 The Crumpling Transition 

The Free energy is now given by 

1 



F(f) = I ci^x 



(13) 



where the dependence on k may be included in the couplings u and v. With 
the leading term having two derivatives, the directions defined by the cou- 
plings u and V are relevant by naive power counting for D < 4. This shows 
that the model is amenable to an e-expansion with e = 4 — D. For practical 
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purposes, it is more convenient to consider the coupling v = v + ^. We pro- 
vide the detailed derivation of the corresponding (3 functions in Appendix 0. 
The result is 

(^u{ur,vr) = -euR + ^{{d/3 + 65/12)u'ji + 6uRVR + A/3vji} , . 
Pv{ur,vr) = -evR + ^{21/16ul + 21/2uRVR + {4d + 5)vl} ^ ^ 

Rather surprisingly, this set of P functions does not possess a FP, except 
for d > 219. This result would suggest that the crumpling transition is first 
order for d = 3. Other estimates, however, give results which are consistent 
with the crumpling transition being continuous. These are 



Limit of large elastic constants ||50||: The Crumpling transition is ap- 



proached from the flat phase, in the limit of infinite elastic constants. 
The model is 



H 



NL 



J d^^a^iArr , (15) 



with the further constraint daTdpr = 5a(3- Remarkably, the /5-function 
may be computed within a large d expansion, yielding a continuous 
crumpling transition with size exponent at the transition (for D = 2) 



SCSA Approximation El : The Schwinger-Dyson equations for the model 
given by Eq.(|13[) are truncated to include up to four point vertices. The 
result for the Hausdorff dimension and size exponent is 



d 



H 



2.73 



0.732 



(17) 



MCRG Calculation |£2[: The crumpling transition is studied using 
MCRG (Monte Carlo Renormalization Group) techniques. Again, the 
transition is found to be continuous with exponents 



d 



H 



2.64(5) 



0.85(9) 



(18) 



Each of these three independent estimates give a continuous crumpling 
transition with a size exponent in the range v ~ 0.7 ± .15. It would be 
interesting to understand how the e-expansion must be performed in order 
to reconcile it with these results. 
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Figure 14: Plot of tlie specific heat observable . The growth of the specific 
heat peak with system size indicates a continuous transition. 



Further evidence for the crumpling transition being continuous is provided 
by numerical simulations ^ 0, [ 



45, 46, 47, 38 1 where the 



analysis of observables hke the specific heat (see Fig.|T4D or the radius of 
gyration radius give textbook continuous phase transitions, although the 
precise value of the exponents at the transition are difficult to pin down. Since 
this model has also been explored numerically with different discretizations 
on several lattices, there is clear evidence for universality of the crumpling 
transition again consistent with a continuous transition. In Appendix P 
we present more details of suitable discretizations of the energy for numerical 
simulations of membranes. 



4.1.3 The Flat Phase 

In a flat membrane (see Fig.|T^), we consider the strain tensor 

Ual3 = daUfj + dfjUa + dahdfjh . 

The free energy Eq.(0) becomes 



F(u, h)= j rf^x 



(19) 



(20) 
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Figure 15: Coordinates for fluctuations in the flat phase 



where we have dropped irrelevant terms. One recognizes the standard Landau 
Free energy of elasticity theory, with Lame coefficients fi and A, plus an 
extrinsic curvature term, with bending rigidity k. These couplings are related 
to the original ones in Eq.(|^ hj ^ = u(^~^, A = 2v('^~^, k = k(^~^ and 
t = -4(/x + f A)C^-2. 

The large distance properties of the flat phase for crystalline membranes 
are completely described by the Free energy Eg. (P(]|) . Since the bending rigid- 
ity may be scaled out at the crumpling transition, the free energy becomes a 
function of ^ and ^. The /3-function for the couplings u, v ai k, = 1 may be 
calculated within an e-expansion, which we describe in detail in Appendix |^. 
Let us recall that the dependence on k may be trivially restored at any stage. 
The result is 



P,{fiR,XR) = -e/i« + |^(| + 20A) (21) 



Pxil^R, Aif) = -eXn + ^{^l^l + 2(4 + 10 A) XRfiR + 2d,Xl) , 

where = d — D, and A = + a ^ • These /3 functions show four fixed points 
whose actual values are shown in Table |I]. 



As apparent from Fig. 16, the phase diagram of the fiat phase turns out 



to be slightly more involved than the one shown in Fig.|l3|, as there are three 
FPs in addition to the FLFP already introduced. These additional FPs are 
infra-red unstable, however, and can only be reached for very specific values 
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Figure 16: Phase diagram for the phantom flat phase. There are three infra- 
red unstable FPs, labelled by FLPl, FLP2 and FLP3, but the physics of the 
flat phase is governed by the infra-red stable FP (FLFP). 

of the Lame coefficients, so for any practical situation we can regard the 
FLFP as the only existing FP in the flat phase. 

The properties of the fiat phase 

The flat phase is a very important phase as will be apparent once we 
study the full model, including self- avoidance. For that reason we turn now 
to a more detailed study of its most important properties. 

Fig.|Tl| (c) gives an intuitive visualization of a crystalline membrane in the 
flat phase. The membrane is essentially a flat two dimensional object up to 
fluctuations in the perpendicular direction. The rotational symmetry of the 
model is spontaneously broken, being reduced from 0{d) to 0{d—D) xO{D). 
The remnant rotational symmetry is realized in Eq.(pO|) as 

/i,(x) ^ /i,(x) + A*"x„ (22) 
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Table 1: The FPs and critical exponents of the flat phase. 



where A^°' is a. D x {d — D) matrix. This relation is very important as it 
provides Ward identities which simplify enormously the renormalization of 
the theory. 

Let us first study the critical exponents of the model. There are two key 
correlators, involving the in-plane and the out-of-plane phonon modes. Using 
the RG equations, it is easy to realize that at any given FP, the low-p limit 
of the model is given by 

Tuuip) ~ bf^"" (23) 

where the last equation defines the anomalous elasticity k{p} as a function 
of momenta p. These two exponents are not independent, since they satisfy 
the scaling relation 

r]u = A-D-2r], (24) 

which follows from the Ward identities associated with the remnant rotational 
symmetry (Eqn.(|2^). Another important exponent is the roughness exponent 
(, which measures the fluctuations transverse to the fiat directions. It can 
be expressed as C = ^— f^. 

The long wavelength properties of the fiat phase are described by the 
FLFP (see Fig.|l^. Since the FLFP occurs at non-zero renormalized values 
of the Lame coefficients, the associated critical exponents discussed earlier 
are clearly non-Gaussian. Within an e-expansion, the values for the critical 
exponents are given in Table There are alternative estimates available 
from different methods. These are 



Numerical Simulation: In ||38[ a large scale simulation of the model 



was performed using very large meshes. The results obtained for the 
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critical exponents are very reliable, namely 



r/„ = 0.50(1) T] = 0.750(5) C = 0.64(2) (25) 
For a review of numerical results see 

SCSA Approximation: This consists of suitably truncating the Schwinger- 



Dyson equations to include up to four-point correlation functions ||5T 
The result for general d is 

^^^^ = 4 + (16-24 + c^g)V^ ' ^''^ 
which for d = 3 gives 

riu = 0.358 T] = 0.821 C = 0.59 (27) 
• Large d expansion: The result is |Q 

V = l^ 7/(3) = 2/3 (28) 

We regard the results of the numerical simulation as our most accurate es- 
timates, since we can estimate the errors. The results obtained from the 
SCSA, which are the best analytical estimate, are in acceptable agreement 
with simulations. 

Finally there are two experimental measurements of critical exponents for 
the flat phase of crystalline membranes. The static structure factor of the 
red blood cell cytoskeleton (see Sect.|ip has been measured by small-angle 
x-ray and light scattering, yielding a roughness exponent of ( = 0.65(10) 



r3| . Freeze- fracture electron microscopy and static light scattering of the 
conformations of graphitic oxide sheets (Sect.p revealed fiat sheets with a 
fractal dimension dn = 2.15(6). Both these values are in good agreement 
with the best analytic and numerical predictions, but the errors are still too 
large to discriminate between different analytic calculations. 

The Poisson ratio of a crystalline membrane (measuring the transverse 
elongation due to a longitudinal stress |5^) is universal and within the SCSA 
approximation, which we regard as the more accurate analytical estimate, is 
given by 

^(D) = --^ ^ a(2) = -1/3 , (29) 
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This result has been accurately checked in numerical simulations |^|. Rather 
remarkably, it turns out to be negative. Materials having a negative Poisson 
ratio are called auxetic. This highlights potential applications of crystalline 
membranes to materials science since auxetic materials have a wide variety 
of potential applications as gaskets, seals etc. 

Finally another critical regime of a flat membrane is achieved by subject- 
ing the membrane to external tension [^6|, This allows a low temperature 
phase in which the membrane has a domain structure, with distinct domains 
corresponding to fiat phases with different bulk orientations. This describes, 
physically, a buckled membrane whose equilibrium shape is no longer planar. 



4.2 Self-avoiding 

Self-avoidance is a necessary interaction in any realistic description of a crys- 
talline membrane. It is introduced in the form of a delta-function repulsion 
in the full model Eq.(|^). We have already analyzed the phantom case and 
explored in detail the distinct phases. The question before us now is the 
effect of self-avoidance on each of these phases. 

The first phase we analyze is the fiat phase. Since self-intersections are 
unlikely in this phase, it is intuitively clear that self-avoidance should be 
irrelevant. This may also be seen if one neglects the effects of the in-plane 
phonons. In the self-avoiding term for the fiat phase we have 

f(y)) (30) 

y) + u(x)-u(y))5'^-^(Mx)-My)) 

y))5^-^(Mx) - My)) = , 

as the trivial contribution where the membrane equals itself is eliminated by 
regularization. The previous argument receives additional support from nu- 
merical simulations in the flat phase, where it is found that self-intersections 
are extremely rare in the typical configurations appearing in those simula- 
tions. It seems clear that self-avoidance is most likely an irrelevant operator, 
in the RG sense, of the FLFP. Nevertheless, it would be very interesting if one 
could provide a more rigorous analytical proof for this statement. A rough 
argument can be made as follows. Shortly we will see that the Flory approx- 
imation for self-avoiding membranes predicts a fractal dimension dn = 2.5. 



-J rf^xdV'(r(x)- 
= ^/rf^xdV^(C(x- 
~ Urf^xrfV^(C(x- 
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For bulk dimension d exceeding 2.5 therefore we expect self-avoidance to be 
irrelevant. A rigorous proof of this sort remains rather elusive, as it involves 
the incorporation of both self-avoidance and non-linear elasticity, and this 
remains a difficult open problem. 

The addition of self-avoidance in the crumpled phase consists of adding 
the self-avoiding interaction to the free energy of Eq. ([T0| ) 

F{r) = y rf^x(9„f(x))2 + y rf^xt/V'(r(x) - f(y)) , (31) 

which becomes the natural generalization of the Edwards' model for polymers 
to D-dimensional objects. Standard power counting shows that the GFP of 
the crumpled phase is infra-red unstable to the self-avoiding perturbation for 

e{D,d) = 2D - d'^-^ > , (32) 

which implies that self-avoidance is a relevant perturbation for D = 2-objects 
at any embedding dimension d. The previous remarks make it apparent that 
it is possible to perform an e-expansion of the model [58, 59, |6^, 31 1. In 
Appendix ^ we present the calculation of the /3-function at lowest order 
in e using the MOPE (Multi-local-operator-product-expansion) formalism 



62| , The MOPE formalism has the advantage that it is more easily 

generalizable to higher orders in e, and enables concrete proofs showing that 
the expansion may be carried out to all orders. At lowest order, the result 
for the /3-function is 



T{^) ^ 2 2D 



hi 
2 



= -ebn + a,bj,. (33) 

The infra-red stable FP is given at lowest order in e by = ^, which 
clearly shows that the GFP of the crumpled phase is infra-red unstable in 
the presence of self-avoidance. 

The preceding results are shown in Fig.|l^ and may be summarized as 

• The flat phase of self-avoiding crystalline membranes is exactly the 
same as the flat phase of phantom crystalline tethered membranes. 

• The crumpled phase of crystalline membranes is destabilized by the 
presence of any amount of self- avoidance. 
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Figure 17: The addition of Self-avoidance at the Crumpled and Flat phases. 

The next issue to elucidate is whether this new SAFP describes a crum- 
pled self-avoiding phase or a flat one and to give a more quantitative descrip- 
tion of the critical exponents describing the universality class. Supposing 
that the SAFP is, in fact, fiat we must understand its relation to the FLFP 
describing the physics of the fiat phase and the putative phase transitions 
between these two. 

4.2.1 The nature of the SAFP 

Let us study in more detail the model described in Eg. (|3T|) . The key issue 
is whether this model still admits a crumpled phase, and if so to determine 
the associated size exponent. On general grounds we expect that there is a 
critical dimension dc-, below which there is no crumpled phase (see Fig.|l^). 

An estimate for the critical dimension may be obtained from a Flory 
approximation in which minimizes the free energy obtained by replacing both 
the elastic and self-avoiding terms with the radius of gyration raised to the 
power of the appropriate scaling dimensions. Within the Flory treatment a 
Z)-dimensional membrane is in a crumpled phase, with a size exponent given 
by 

v={D + 2)/{d + 2) . (34) 

From this it follows that dc = D (see Fig.|l8D. The Flory approximation, 
though very accurate for polymers {D = 1), remains an uncontrolled ap- 
proximation. 

In contrast the e-expansion provides a systematic determination of the 
critical exponents. For the case of membranes, however, some extrapolation 
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Figure 18: The size exponent as a function of d. There is a critical dimension 
dc below which the crumpled phase disappears. 



is required, as the upper critical dimension is infinite. This was done in |]6J] , 
where it is shown that reconsidering the e-expansion as a double expansion 
in e and -D, critical quantities may be extrapolated for D = 2-dimensional 
objects. At lowest order in e, the membrane is in a crumpled phase. The 
enormous task of calculating the next correction (e^) was successfully carried 
out in employing more elaborate extrapolation methods than those in 
|6^ . Within this calculation, the d = 3 membrane is still in a crumpled 
phase, but with a size exponent now closer to 1. It cannot be ruled out that 
the e-expansion, successfully carried out to all orders could give a flat phase 
V = 1. In fact, the authors in [^, ^] present some arguments in favor of a 
scenario of this type, with a critical dimension dc ~ 4. 

Other approaches have been developed with different results. A Gaussian 
approximation was developed in p^, |63] . The size exponent of a self-avoiding 



(35) 



membrane within this approach is 



4/d, 



and since one has v > 1 for (i < 4, one may conclude that the membrane 
is fiat for c? < = 4. Since we cannot determine the accuracy of the 
Gaussian approximation this estimate must be viewed largely as interesting 
speculation. Slightly more elaborate arguments of this type yield an 
estimated critical dimension d^ = 3. 



29 



Numerical simulations 



We have seen that numerical simulations provide good support for an- 
alytic results in the case of phantom membranes. When self-avoidance is 
included, numerical simulations become invaluable, since analytic results are 
harder to come by. It is for this reason that we discuss them in greater detail 
than in previous sections. 

A possible discretization of membranes with excluded volume effects con- 
sists of a network of N particles arranged in a triangular array. Nearest 
neighbors interact with a potential 

T/ Jo for \ f\<b , . 

= I oo for |r1 > 6 ' ^^^^ 

although some authors prefer a smoothened version, with the same general 
features. The quantity b is of the order of a few lattice spacings. This is 
a lattice version of the elastic term in Eg. (pi]) . The discretization of the 
self-avoidance is introduced as a repulsive hard sphere potential, now acting 
between any two atoms in the membrane, instead of only nearest neighbors. 
A hard sphere repulsive potential is, for example, 

where a is the range of the potential, and a < b. Again, some smoothened 
versions, continuous at |r| = a, have also been considered. This model may 
be pictured as springs, defined by the nearest-neighbor potential Eg. (p6D , 
with excluded volume effects enforced by balls of radius a (Eg. (|37D ). This 
model represents a lattice discretization of Eg. . 

Early simulations of this type of model [^9|, ^ provided a first estimate of 
the size exponent at d = 3 fully compatible with the Flory estimate Eg. (^4]) . 
The lattices examined were not very large, however, and subseguent simula- 
tions with larger volumes pUl I7D| found that the d = 3 membrane is actually 



fiat. This result is even more remarkable if one recalls that there is no explicit 
bending rigidity. 

The fiat phase was a very surprising result, in some conflict with the 
insight provided from the analytical estimates discussed in the previous sub- 



section. An explanation for it came from the observation that excluded 



volume effects induce bending rigidity, as depicted in Fig. 19. The reason is 
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Figure 19: Visualization of bending rigidity generated by a hard sphere po- 
tential Normals ni and n2 cannot be anti-parallel. 



that the excluded volume effects generate a non-zero expectation value for 
the bending rigidity, since the normals can be parallel, but not anti-parallel 



see 



191) . This induced bending rigidity was estimated and found to be big 



enough to drive the self-avoiding membrane well within the flat phase of the 
phantom one. This means that this particular discretization of the model 
renders any potential SAFP inaccessible and the physics is described by the 
FLFP. In ||72[ the structure function of the self-avoiding model is numerically 
computed and found to compare well with the analytical structure function 
for the fiat phase of phantom crystalline membranes, including comparable 
roughness exponents. 

The natural question then to ask is whether it is possible to reduce the 
bending rigidity sufficiently to produce a crumpled self-avoiding phase. Sub- 
sequent studies addressed this issue in various ways. The most natural way 
is obviously to reduce the range of the potential sufficiently that the induced 
bending rigidity is within the crumpled phase. This is the approach followed 



m 



73|. The flat phase was found to persist to very small values of a, with 
eventual signs of a crumpled phase. This crumpled phase may essentially 
be due to the elimination at self-avoidance at sufficiently small cr. A more 
comprehensive study, in which the same limit is performed this time with 
an excluded volume potential which is a function of the internal distance 
along the lattice ||7^, concluded that for large membranes, inclusion of ex- 
cluded volume effects, no matter how small, leads to flatness. A different 
approach to weakening the flat phase, bond dilution ||75|, found that the flat 
phase persists until the percolation critical point. In conclusion the bulk 
of accumulated evidence indicates that flatness is an intrinsic consequence 
of self-avoidance. If this is indeed correct the SAFP coincides with FLFP 
and this feature is an inherent consequence of self-avoidance, rather than an 
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artifact of discretization. 

Given the difficulties of finding a crumpled phase with a repulsive poten- 
tial, simulations for larger values of the embedding space dimension d have 



also been performed ^ . These simulations show clear evidence that the 
membrane remains fiat for d = 3 and 4 and undergoes a crumpling transition 
for (i > 5, implying dc> 4:. 

An alternative approach to incorporating excluded volume effects cor- 
responds to discretize a surface with a triangular lattice and imposing the 
self-avoidance constraint by preventing the triangular plaquettes from inter- 
penetrating. This model has the advantage that is extremely fiexible, since 
there is no restriction on the bending angle of adjacent plaquettes (triangles) 
and therefore no induced bending rigidity (see y). 

The first simulations of the plaquette model found a size exponent 
in agreement with the Flory estimate Eq. A subsequent simulation f^^ 



disproved this result, and found a size exponent u = 0.87, higher than the 
Flory estimate, but below one. More recent results using larger lattices 
and more sophisticated algorithms seem to agree completely with the results 
obtained from the ball and spring models |^ . 

Further insight into the lack of a crumpled phase for self-avoiding crys- 



talline membranes is found in the study of folding ||8T|, q2|, |8^, ^ |85|, ^y, |87 
This corresponds to the limit of infinite elastic constants studied by David 
and Guitter with the further approximation that the space of bending angles 
is discretized. One quickly discovers that the refiection symmetries of the 
allowed folding vertices forbid local folding (crumpling) of surfaces. There is 
therefore essentially no entropy for crumpling. There is, however, local un- 
folding and the resulting statistical mechanical models are non-trivial. The 
lack of local folding is the discrete equivalent of the long-range curvature- 
curvature interactions that stabilize the fiat phase. The dual effect of the 
integrity of the surface (time-independent connectivity) and self- avoidance is 
so powerful that crumpling seems to be impossible in low embedding dimen- 
sions. 

4.2.2 Attractive potentials 

Self- avoidance, as introduced in Eq. (p7|) is a totally repulsive force among 
monomers. There is the interesting possibility of allowing for attractive po- 
tentials also. This was pioneered in |7l|] as a way to escape to the induced 



bending rigidity argument (see Fig.O), since an attractive potential would 
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correspond to a negligible (or rather a negative) bending rigidity. Remark- 
ably, in , a compact (more crumpled) self-avoiding phase was found, with 



fractal dimension close to 3. 

This was further studied in [Q, where it was found that with an attractive 
Van der Waals potential, the crystalline membrane underwent a sequence of 
folding transitions leading to a crumpled phase. In ^| similar results were 



found, but instead of a sequence of folding transitions a crumpled phase 
was found with an additional compact (more crumpled phase) at even lower 



temperatures. Subsequent work gave some support to this scenario ||90 



On the analytical side, the nature of the G-point for membranes and its 



relevance to the issue of attractive interactions has been addressed in ^1 . 

We think that the study of a tether with an attractive potential remains 
an open question begging for new insights. A thorough understanding of the 
nature of the compact phases produced by attractive interactions would be 
of great value. 

4.2.3 The properties of the SAFP 

The enormous efforts dedicated to study the SAFP have not resulted in 
a complete clarification of the overall scenario since the existing analytical 
tools do not provide a clear picture. Numerical results clearly provide the best 
insight. For the physically relevant case d = 3, the most plausible situation is 
that there is no crumpled phase and that the flat phase is identical to the fiat 
phase of the phantom model. For example, the roughness exponents (sa from 
numerical simulations of self-avoidance at d = 3 using ball-and-spring models 
76| and the roughness exponent at the FLFF, Eg . (^5]) , compare extremely 



well 

CsA = 0.64(4) , C = 0.64(2) , (38) 

So the numerical evidence allows us to conjecture that the SAFP is exactly 
the same as the FLFP, and that the crumpled self-avoiding phase is absent in 
the presence of purely repulsive potentials (see Figj20|). This identification of 
fixed points enhances the significance of the FLFP treated earlier. It would 
be very helpful if analytical tools were developed to further substantiate this 
statement. 
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Figure 20: The conjectured phase diagram for self-avoiding crystalhne mem- 
branes in d = 3. With any degree of self-avoidance the flows are to the flat 
phase fixed point of the phantom model (FLFP). 

5 Anisotropic Membranes 

An anisotropic membrane is a crystalline membrane having the property that 
the elastic or the bending rigidity properties in one distinguished direction 
are different from those in the D — 1 remaining directions. As for the isotropic 
case we keep the discussion general and describe the membrane by a d- 
dimensional r(xx,|/), where now the D dimensional coordinates are split 
into D — 1 coordinates and the orthogonal distinguished direction y. 

The construction of the Landau free energy follows the same steps as in 
the isotropic case. Imposing translational invariance, 0{d) rotations in the 
embedding space and 0{D — 1) rotations in internal space, the equivalent of 
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Eq.(^ is now 



F(f(x)) 



+ f^±yd'yf- dlf+ + tyidyf)^ 



+v±y{d^rr{dyrr 




(39) 



This model has eleven parameters, representing distinct physical interactions: 

• i^±,i^y,i^±y bending rigidity: the anisotropic versions of the isotropic 
bending rigidity splits into three distinct terms. 

• t±,ty,u±±,Uyy,v±±,v±y clastlc constants : there are six quantities de- 
scribing the microscopic elastic properties of the anisotropic membrane. 

• b, self-avoidance coupling: This particular term is identical to its isotropic 
counterpart. 

Following the same steps as in the isotropic case, we split 



with u_L being the D — 1-dimensional phonon in-plane modes, Uy the in-plane 
phonon mode in the distinguished direction y and h the d — D out-of-plane 
fluctuations. If C± = Cy = 0, the membrane is in a crumpled phase and if 
both C± 7^ and Cy 7^ the membrane is in a fiat phase very similar to the 
isotropic case (how similar will be discussed shortly). There is, however, the 
possibility that C± = and Cy 7^ or C± 7^ and (y = 0. This describes a 
completely new phase, in which the membrane is crumpled in some internal 
directions but fiat in the remaining ones. A phase of this type is called a 
tubular phase and does not appear when studying isotropic membranes. In 
Fig.^ we show an intuitive visualization of a tubular phase along with the 
corresponding fiat and crumpled phases of anisotropic membranes. 

We will start by studying the phantom case. We show, using both analj^- 
ical and numerical arguments, that the phase diagram contains a crumpled. 



r(x) = (C±Xi + Ui(x), (yy + Uy{x.), h{x.)) 



(40) 
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Figure 21: Examples of a) the tubular phase b) the crumpled phase and c) 
the flat phase of an anisotropic phantom crystalline membrane taken from 
the simulations of |^| 



tubular and flat phase. The crumpled and flat phases are equivalent to the 
isotropic ones, so anisotropy turns out to be an irrelevant interaction in those 
phases. The new physics is contained in the tubular phase, which we describe 
in detail, both with and without self-avoidance. 

5.1 Phantom 

5.1.1 The Phase diagram 

We first describe the mean field theory phase diagram and then the effect 
of fluctuations. There are two situations depending on the particular values 
of the function A, which depends on the elastic constants u±±, v±y, Uyy and 
f_L_L. Since the derivation is rather technical, we refer to Appendix ^ for the 
details. 

• Case A (A > 0): the mean field solution displays all possible phases. 
When ty > and t± > the model is in a crumpled phase. Lowering 
the temperature, one of the t couplings becomes negative, and we reach 
a tubular phase (either ± or y-tubule). A further reduction of the 
temperature eventually leads to a flat phase. 
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Case B (A < 0): in this case the flat phase disappears from the mean 
field solution. Lowering further the temperature leads to a continuous 
transition from the crumpled phase to a tubular phase. Tubular phases 
are the low temperature stable phases in this regime. 



This mean field result is summarized in Fig.|22 



A) A>0 



B) A<0 



CRUMPLED 
PHASE 




CRUMPLED 
PHASE 




Figure 22: The phase diagram for anisotropic phantom membranes 

Beyond mean field theory, the Ginsburg criterion applied to this model 
tells us that the phase diagram should be stable for physical membranes 
D = 2 at any embedding dimension d, so the mean field scenario should give 
the right qualitative picture for the full model. 

Numerical simulations have spectacularly confirmed this result. We have 
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already shown in Fig.^ the results from the numerical simulation in 
where it was shown that changing the temperature generates a sequence of 
transitions crumpled-to-tubular and tubular-to-flat, in total agreement with 



case A) in the mean field result illustrated in Fig.^. 

We now turn to a more detailed study of both the crumpled and flat 
anisotropic phases. Since we have already studied crumpled and flat phases 
we just outline how those are modified when anisotropy is introduced. 
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5.1.2 The Crumpled Anisotropic Phase 

In this phase ty > and t± > 0, and the free energy Eg. (|39|) reduces for 
> 2 to 

F(f(x)) = ^J d^'Wdy [t±{d^r)^ + ty{dyr)^] + Irrelevant . (41) 

By redefining the y direction as y' = j^y this reduces to Eg. (p!OD , with t = t±. 
We have proved that anisotropy is totally irrelevant in this particular phase. 

5.1.3 The Flat Phase 

This phase becomes eguivalent to the isotropic well. Intuitively, this 

may be obtained from the fact that if the membrane is flat, the intrinsic 
anisotropics are only apparent at short-distances, and therefore by analyzing 
the RG flow at larger and larger distances the membrane should become 



isotropic. This argument may be made slightly more precise ||9^ . 



5.2 The Tubular Phase 

We now turn to the study of the novel tubular phase, both in the phantom 
case and with self-avoidance. Since the physically relevant case for mem- 
branes is Z) = 2 the y-tubular and ±-tubular phase are the same. So we 
concentrate on the properties of the y-tubular phase. 

The key critical exponents characterizing the tubular phase are the size 
(or Flory) exponent z/, giving the scaling of the tubular diameter Rg with the 
extended {Ly) and transverse {L±) sizes of the membrane, and the rough- 
ness exponent ( associated with the growth of height fluctuations hrms (see 
Fig.pp: 



Rg{L^,Ly) oc L'iSniLy/Ll) (42) 
hrmsiL±, Ly) (X LySh{Ly/Lj_) 



Here Sr and Sh are scaling functions @, |95| and z is the anisotropy exponent. 

The general free energy described in Eg.(^) may be simplifled consider- 
ably in a y-tubular phase. The analysis reguired is involved and we refer the 
interested reader to 0. We just guote the flnal result. It is 



F{u, h) = d^~^x±dy [Kid^hf + t{dji) 
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Figure 23: A schematic illustration of a tubular configuration indicating the 
radius of gyration Rg and the height fluctuations hrms- 

+g±{daU + dahdyhf 
+ 9y{dyU + \{dyhf? 

+^ j dyd^-'^^d^-'^j'-'iK^^^y) - hi^'^,y)) . (43) 

Comparing with Eq.(^), this free energy does represent a simplification as 
the number of couplings has been reduced from eleven to five. Furthermore, 
the coupling g_i_ is irrelevant by standard power counting. The most natural 
assumption is to set it to zero. In that case the phase diagram one ob- 
tains is shown in Fig.^. Without self-avoidance 6 = 0, the Gaussian Fixed 
Point (GFP) is unstable and the infra-red stable FP is the tubular phase FP 
(TPFP). Any amount of self-avoidance, however, leads to a new FP, the Self- 
avoiding Tubular FP (SAFP), which describes the large distance properties 
of self-avoiding tubules. 

We just mention, though, that other authors advocate a different scenario 



95| . For sufficiently small embedding dimensions d, including the physical 
d = 3 case, these authors suggest the existence of a new bending rigidity 
renormalized FP (BRFP), which is the infra-red FP describing the actual 
properties of self-avoiding tubules (see Fig. B^). 



Here we follow the arguments presented in ||9^ and consider the model 
defined by Eq. (l43|) with the g±-teTm as the model describing the large dis- 
tance properties of tubules. One can prove then than there are some general 
scaling relations among the critical exponents. All three exponents may be 
expressed in terms of a single exponent 

3 1-D 

c 



2 2z 

= (z . (44) 
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Figure 24: The phase diagram for self-avoiding anisotropic membranes with 
the Gaussian fixed point (GFP), the tubular phase fixed point (TPFP) and 
the self-avoidance fixed point (SAFP). 



Remarkably, the phantom described by Eq.(H3) can be solved 

exactly. The result for the size exponent is 

^Phantom ) ^Phantom (2) = 4 (45) 



with the remaining exponents following from the scaling relations Eq.( 

The self-avoiding case may be treated with techniques similar to those in 
isotropic case. The size exponent may be estimated within a Flory approach. 
The result is 

D + l , , 

^Flory = -J-^ ■ (46) 

The Flory estimate is an uncontrolled approximation. Fortunately, a e- 
expansion, adapting the MOPE technique described for the self-avoiding 
isotropic case to the case of tubules, is also possible 0, ^ . The /^-functions 



are computed and provide evidence for the phase diagram shown in Fig.|2^. 
Using rather involved extrapolation techniques, it is possible to obtain esti- 



40 




Figure 25: The phase diagram for self-avoiding anisotropic membranes with 
the Gaussian fixed point (GFP), the tubular phase fixed point (TPFP), 
the self-avoidance fixed point (SAFP) and the bending rigidity fixed point 
(BRFP). 

mates for the size exponent, which are shown in Table ^. The rest of the 
exponents may be computed from the scaling relations. 

6 Defects in membranes: The Crystalline- 
Fluid transition and Fluid membranes 

A flat crystal melts into a liquid when the temperature is increased. This 
transition may be driven by the sequential liberation of defects, as predicted 
by the KTNHY theory. The KTNHY theory is schematically shown in Fig.|26|. 
With increasing temperature, a crystal melts first to an intermediate hexatic 
phase via a continuous transition, and finally goes to a conventional isotropic 
fluid phase via another continuous transition. 

We will not review here either the KTNHY theory or the experimental 



evidence in its favor - 98 1. We just want to emphasize here that the KTNHY 



theory is in general agreement with existing experiments, although there are 



41 



d 
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vf 
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^Flory 


8 


0.333(5) 


0.34(1) 


0.34(1) 


0.333 


7 


0.374(8) 


0.39(2) 


0.39(2) 


0.375 


6 


0.42(1) 


0.44(2) 


0.44(4) 


0.429 


5 


0.47(1) 


0.51(3) 


0.51(5) 


0.500 


4 


0.54(2) 


0.60(4) 


0.60(6) 


0.600 


3 


0.62(2) 


0.71(6) 


0.70(9) 


0.750 



Table 2: Different extrapolations for the size exponent at different embedding 
dimensions d The first column gives the corrections to the mean field 



result, the second corrections to the Flory estimate and the third corresponds 
to corrections to the Gaussian approximation. The last column quotes the 
Flory estimate for comparison. 
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Figure 26: Two stage melting according to KTNHY theory. 



two main points worth keeping in mind when studying the more difficult case 
of fluctuating geometries. 1) The experimental evidence for the existence of 
the hexatic phase is not completely settled in those transitions which are 
continuous. 2) Some 2D crystals (like Xenon absorbed on graphite) melt to 
a fluid phase via a first order transition without any intermediate hexatic 
phase. 

The straight-forward translation of the previous results to the tethered 
membrane would suggest a similar scenario. There would then be a crys- 
talline to hexatic transition and a hexatic to fluid transition, as schemati- 
cally depicted in Fig.|^. Although the previous scenario is plausible, there 
are no solid experimental or theoretical results that establish it. From the 
theoretical point of view, for example, an important open problem is how to 
generalize the RG equations of the KTNHY theory to the case of fluctuating 
geometry. The situation looks even more uncertain experimentally, especially 
considering the elusive nature of the hexatic phase even in the case of flat 
monolayers. 

In this review we will assume the general validity of the KTNHY scenario 
and we describe models of hexatic membranes, as well as fluid membranes. 
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The study of the KTNHY theory in fluctuating geometries is a fascinating 
and challenging problem that deserves considerable effort. In this context, let 
us mention recent calculations of defects on frozen topographies |^^, which 
show that even in the more simplified case when the geometry is frozen, de- 
fects proliferate in an attempt to screen out Gaussian curvature, even at zero 
temperature, and organize themselves in rather surprising and unexpected 
structures. These results hint at a rich set of possibilities for the more general 
case of fluctuating geometries. 

6.1 Topological Defects 

A crystal may have different distortions from its ground state. Thermal 
fluctuations are the simplest. Thermal fluctuations are small displacements 
from the ground state, and therefore one may bring back the system to its 
original positions by local moves without affecting the rest of the lattice. 
There are more subtle lattice distortions though, where the lattice cannot be 
taken to its ground state by local moves. These are the topological defects. 
There are different possible topological defects that may occur on a lattice, 
but we just need to consider dislocations and disclinations. Let us review 
the most salient features. 




Figure 27: Example of a dislocation showing the breaking of the translational 
holonomy measured by the Burgers vector b. 

• Dislocation : represents the breaking of the translational holonomy. A 
path that would naturally close in a perfect lattice fails to close by a 
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vector b, the Burgers vector, as illustrated in Fig.^. In a flat mono- 
layer, the energy is 



E 



Sir ^a' 



(47) 



where Kq is the Young Modulus. It diverges logarithmically with system 
size. 

Disclination : represents the breaking of the rotational holonomy. The 
bond angle around the point defect is a multiple of the natural bond 
angle in the ground state (| in a triangular lattice), as illustrated in 
Fig.^ for a + and — disclination. The energy for a disclination in a 
flat monolayer is given by 



E 



327r 



(48) 



Note the quadratic divergence of the energy with system size R. 




Figure 28: Example of a minus disclination (left figure) and a plus disclination 
(right figure). The orientational holonomy is broken by ±| respectively. 

Inspection of Fig.^ shows that a dislocation may be regarded as a tightly 
bound +,- disclination pair. 

6.1.1 Topological Defects in fluctuating geometries 

The problem of understanding topological defects when the geometry is al- 
lowed to fluctuate was addressed in ||10CI|| (see |p.01|| for a review). The impor- 
tant new feature is that the energy of a disclination defect may be lowered 
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considerably if the membrane buckles out-of-the-plane. That is, the mem- 
brane trades elastic energy for bending rigidity. The energy for a buckled 
free disclination is given by 

E = f{K,K,,q,)H-) , (49) 

where / is some complicated function that may be evaluated numerically for 
given values of the parameters. It depends explicitly on g^, which implies 
that the energies for positive and negative disclinations may be different, 
unlike the situation in flat space. The extraordinary reduction in energy 
from to InR is possible because the buckled membrane creates positive 
Gaussian curvature for the plus-disclination and negative curvature for the 
negative-disclination. This is a very important physical feature of defects 
on curved surfaces. The defects attempt to screen out like-sign curvature, 
and analogously, like-sign defects may force the surface to create like-sign 
curvature in order to minimize the energy. 

The reduction in energy for a dislocation defect is even more remarkable, 
since the energy of a dislocation becomes a constant, independent of the 
system size, provided the system is larger than a critical radius Re- Again, 
by allowing the possibihty of out-of-plane buckling, a spectacular reduction 
in energy is achieved (from R^ to In R) . 

The study of other topological defects, e.g. vacancies, interstitials, and 
grain boundaries, may be carried out along the same lines. Since we are not 



going to make use of it, we refer the reader to the excellent review in ||101j|. 



6.1.2 Melting and the hexatic phase 

The celebrated Kosterlitz-Thouless argument shows that defects will neces- 
sarily drive a 2D crystal to melt. The entropy of a dislocation grows loga- 
rithmically with the system size, so for sufficiently high temperature, entropy 
will dominate over the dislocation energy (Eq.(|47|)) and the crystal will nec- 
essarily melt. If the same Kosterlitz-Thouless is applied now to a tethered 
membrane, the entropy is still growing logarithmically with the system size, 
while the energy becomes independent of the system size, as explained in the 
previous subsection, so any finite temperature drive the crystal to melt, and 
the low temperature phase of a tethered membrane will necessarily be a fiuid 
phase, either hexatic if the KTNHY melting can be applied, or a conventional 
fiuid if a first order transition takes place, or even some other more perverse 
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possibility. This problem has also been investigated in numerical simulations 
|102|| , which provide some concrete evidence in favor of the KTNHY scenario, 



although the issue is far from being settled. 

It is apparent from these arguments that a hexatic membrane is a very 
interesting and possibly experimentally relevant membrane to understand. 

6.2 The Hexatic membrane 

The hexatic membrane is a fluid membrane that, in contrast to a conven- 
tional fluid, preserves the orientational order of the original lattice (six-fold 
(hexatic) for a triangular lattice). The mathematical description of a fluid 
membrane is very different from those with crystalline order. Since the de- 
scription cannot depend on internal degrees of freedom, the free energy must 
be invariant under reparametrizations of the internal coordinates (that is, 
should depend only on geometrical quantities, or in more mathematical ter- 
minology, must be diffeomorphism invariant). The corresponding free energy 
was proposed by Helfrich | |1U3|| and it is given by 

^ = /i/v^+^/rfxy^^^ (50) 

where n is the bare string tension, k the bending rigidity, g is the determinant 
of the metric of the surface 

5(^^(x) = d^f{x.)d^f{x.) (51) 

and H is the mean Gaussian curvature of the surface. For a good description 
of the differential geometry relevant to the study of membranes we refer 



to ||104|| A hexatic membrane has an additional degree of freedom, the bond 
angle, which is introduced as a field on the surface 6. The hexatic free energy 
||105|| is obtained from adjoining to the fluid case of Eg. (|50|) , the additional 
energy of the bond angle 

nnejT =^J d^^g^'^id^e + n^^ng - ^i){d.O + n^^ng - (52) 

where Ka is called the hexatic stiffness, and is the connection two form 
of the metric, which may be related to the Gaussian curvature of the surface 
by 

K(x) = -^e'^'d^n, (53) 
V9 
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and ^sing is similarly related to the topological defect density 

vr 1 ^ 

s(x) = -— 5]g.5(x,Xi) . (54) 

V y i=i 

From general theorems on differential geometry one has the relation 

i 

V^s(x)=47r-.^g, = 2x, (55) 

where x is the Euler characteristic of the surface. 

Therefore the total free energy for the hexatic membrane is given by 

1^ I" • —'72 



■H/T = f, I I d^^H' + (56) 



2 



The partition function is therefore 



Af+ AT- 

^W(r(x),0(x))/T 



/1=1 u = l 

where y is the fugacity of the disclination density. The partition function 
includes a discrete sum over allowed topological defects, those satisfying the 
topological constraint Eg . (p5D , and a path integral over embeddings r and 
bond angles 6. The previous model remains quite intractable since the sum 
over defects interaction is very difficult to deal with. 

Fortunately, the limit of very low fugacity y —>■ is analytically tractable 
as was shown in the beautiful paper [|107|| . The RG functions can be computed 
within a combined large d and large bending rigidity expansion. The f3 
functions in that limit is given by 

/?(«) = ^ {~Da' + ^a' + 0(1/Kl)) , 
PIKa) = (58) 
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Figure 29: Phase diagram according to ||106|1 . The hexatic membrane inter- 
polates between a crystalhne membrane and a fluid one. For small rigidity 
and large hexatic stiffness the RG flows towards a fluid crumpled phase. For 
small hexatic stiffness and large k it flows towards a tethered like phase, 
whereas for both large k and Ka the flow is to the crinkled phases described 
by the non-trivial FP of Eq. . 



where a = I/k. The physics of hexatic membranes in the limit of very low 
fugacity is very rich and show a line of fixed points parametrized by the 
hexatic stiffness K^- The normal- normal correlation function, for example, 



reads [[TO^] 



(n(r)n(O)) ~ |x|-'' , (59) 

with ?7 = -^d{d — 2)^1^. The FPs of Eg . (|58|) describe a new crinkled phase, 
more rigid than a crumpled phase but more crumpled than a flat one. The 
Hausdorff dimension at the crinkled phase is given by ||107|| 



^ , d{d - 2) ksT 

From the RG point of view, the properties of these crinkled phases are really 
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interesting, since they involve a line of Fixed Points which are inequivalent 
in the sense that the associated critical exponents depend continuously on 
Kai a situation reminiscent of the XF-model. In ||106|| the phase diagram is 
discussed, and the authors propose the scenario depicted in Fig.|29|. How these 
scenarios are modified when the fugacity is considered is not well established 
and we refer the reader to the original papers 



T08i [T09i |T06|1 



The shape fluctuations of hexatic vesicles, for large defect core energies, 



have also been investigated [110, 111 



7 The Fluid Phase 

The study of fluid membranes is a broad subject, currently under intense 
experimental and theoretical work. The Hamiltonian is given by 

n/T = ixl^+^J d^^H' + ^I^K , (61) 

This corresponds to the Helfrich hamiltonian together with a term that allows 
for topology changing interactions. For fixed topology it is well-known [[112| , 
113| , |114| , |115| ] that the one loop beta function for the inverse-bending rigidity 
has a fixed point only at k = 0, which corresponds to the bending rigidity 
being irrelevant at large length scales. The RG flow of the bending rigidity 
is given by 

= Ko- — ln(//a) , (62) 

where a is a microscopic cutoff length. The fluid membrane is therefore crum- 
pled, for arbitrary microscopic bending rigidity kq, at length scales beyond a 
persistence length which grows exponentially with kq. For a fluid membrane 
out-of-plane fluctuations cost no elastic energy (the membrane flows inter- 
nally to accommodate the deformation) and the bending rigidity is therefore 
softened by thermal undulations at all length scales, rather than stiffening 
at long length scales as in the crystalline membrane. 

So far we have assumed an infinite membrane, which is not always a real- 
istic assumption. A thickness may be taken into account via a spontaneous 
extrinsic curvature Hq. The model described by Eq. ^ gets replaced then 
by 

n/T = fiJ^+^J d^^{H -H,Y + ^J^K . (63) 
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Further effects of a finite membrane size for spherical topology have been 
discussed in ||116| , |117| , |118|| . The phase diagram of fluid membranes when 
topology change is allowed is fascinating and not completely understood. A 
complete description of these phases goes beyond the scope of this review - 
we refer the reader to 111, |2^, 



and references therein. 



8 Conclusions 

In this review we have described the distinct universality classes of mem- 
branes with particular emphasis on crystalline membranes. In each case we 
discussed and summarized the key models describing the interactions of the 
relevant large distance degrees of freedom (at the micron scale). The body 
of the review emphasizes qualitative and descriptive aspects of the physics 
with technical details presented in extensive appendices. We hope that the 
concreteness of these calculations gives a complete picture of how to extract 
relevant physical information from these membrane models. 

We have also shown that the phase diagram of the phantom crystalline 
membrane class is theoretically very well understood both by analytical and 
numerical treatments. To complete the picture it would be extremely valu- 
able to find experimental realizations for this particular system. An exciting 
possibility is a system of cross-linked DNA chains together with restriction 
enzymes that catalyze cutting and rejoining The difficult chemistry 



involved in these experiments is not yet under control, but we hope that 
these technical problems will be overcome in the near future. 

There are several experimental realizations of self-avoiding polymerized 
membranes discussed in the text. The experimental results compare very 
well with the theoretical estimates from numerical simulations. As a future 
theoretical challenge, analytical tools need to be sharpened since they fail to 
provide a clear and unified picture of the phase diagram. On the experimental 
side, there are promising experimental realizations of tethered membranes 
which will allow more precise results than those presently available. Among 
them there is the possibility of very well controlled synthesis of DNA networks 
to form physical realizations of tethered membranes. 

The case of anisotropic polymerized membranes has also been described in 
some detail. The phase diagram contains a new tubular phase which may be 
realized in nature. There is some controversy about the precise phase diagram 
of the model, but definite predictions for the critical exponents and other 
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quantities exist. Anisotropic membranes are also experimentally relevant. 
They may be created in the laboratory by polymerizing a fluid membrane in 
the presence of an external electric field. 

Probably the most challenging problem, both theoretically and experi- 
mental, is a complete study of the role of defects in polymerized membranes. 
There are a large number of unanswered questions, which include the exis- 
tence of hexatic phases, the properties of defects on curved surfaces and its 
relevance to the possible existence of more complex phases. This problem is 
now under intense experimental investigation. In this context, let us men- 
tion very recent experiments on Langmuir films in a presumed hexatic phase 



120|| . The coalescence of air bubbles with the film exhibit several puzzling 
features which are strongly related to the curvature of the bubble. 

Crystalline membranes also provide important insight into the fluid case, 
since any crystalline membrane eventually becomes fluid at high temperature. 
The physics of fluid membranes is a complex and fascinating subject in itself 
which goes beyond the scope of this review. We highlighted some relevant 
experimental realizations and gave a quick overview of the existing theoretical 
models. Due to its relevance in many physical and biological systems and its 
potential applications in material science, the experimental and theoretical 
understanding of fluid membranes is, and will continue to be, one of the most 
active areas in soft condensed matter physics. 

We have not been able in this review, simply for lack of time, to address 
the important topic of the role of disorder. We hope to cover this in a separate 
article. 

We hope that this review will be useful for physicists trying to get a 
thorough understanding of the fascinating field of membranes. We think it 
is a subject with significant prospects for new and exciting developments. 

Note: The interested reader may also find additional material in a forthcom- 



ing review by Wiese [|121||, of which we have seen only the table of contents. 
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A Useful integrals in dimensional regulariza- 
tion 

In performing the e-expansion, we will be considering integrals of the form 
/.......(a,&)(p1=/c^"gf^^^, (64) 

These integrals may be computed exactly for general a, h and a = 1, ■ ■ ■ , A^. 
The result will be published elsewhere. We will content ourselves by quoting 
what we need, the poles in e, for the integrals that appear in the diagram- 
matic calculations. We just quote the results 

-^aia2 (2) 2) = — „ APctiPa2~ (65) 
-^010203(2)2) = g^2pA P<^iP<^2Pa3 ~ (66) 
-^«i«2a3a4 (2) 2) — - J {Pa\Pa2Pa3P<^i 

\ 1 

TTT (^<^ic*2 '^0304 ~l~ ^ a\a3^ C12CIA ~l~ '^QiQ4'^a2CK3 ) I (6"^) 

24 e 



-^010203(2; 1) o 2 2 ( (Pcn'^0203 Pa2^a\az ~l~Po3^a2CH 



1 

-PaiPaaPaa)- (68) 



e 



B Some practical identities for RG quantities 



The beta functions defined in Eq. ^ may be re-expressed as 
The previous expression may be further simplified noticing 



(69) 



9 In M 9 In u 
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so that 
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t^i? I 
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(71) 



D 



+ 



where the last result follows from Taylor-expanding. These formulas easily 
allow to compute the corresponding /3-functions. If 



u 



Ur + -(aiiM^ + auURVR + aisw^) 

Vr + -(a2lM^ + a22URVR + a23f^) 



(72) 



from Eq. |7T| and Eq. ^ we easily derive the leading two orders in the couplings 

Pu{ur,vr) = -euR + anu% + ai2UR,VR + aisvj^ (73) 

Pv{ur, Vr) = -EVr + a2iu\ + a22URVR + a23^fi • 

The formula for 7 in Eq. ^ may also be given a more practical expression. 
It is given by 

l={(^uj:— + ^j:—)\riZ^, (74) 
OUR dvR 

Those are the formulas we need in the calculations we present in this 
review. 



C Discretized Model for tethered membranes 

In this appendix we present appropriate discretized models for numerical 
simulation of tethered membranes. The surface is discretized by a triangu- 
lar lattice defined by its vertices {r}a=i,-, with a corresponding discretized 
version of the Landau elastic term Eq. ^ given by | |100| | 

Fs = ^Y.i^\ra-n\~l)\ (75) 

where (a, h) are nearest-neighbor vertices. If we write ra = + Uq with 
defining the vertices of a perfectly regular triangular lattice and u the small 
perturbations around it, one gets 

\ra -n\ = l + + • • • , (76) 
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with Uap being a discretized strain tensor and we we have neglected higher 
order terms. Plugging the previous expression into Eq. |75| and passing from 
the discrete to the continuum language we obtain 



13 I £-M2u%, + ulJ (77) 



8 

which is the elastic part of the free energy Eq. 20 with X = fi = ^(3. 
The bending rigidity term is written in the continuum as 

Se,t = J (f^^KpKf = J £vi^g^''Vc.nVpn (78) 

where n is the normal to the surface and V is the covariant derivative (see 



122|| for a detailed description of these geometrical quantities). We discretize 



the normals form the the previous equation by 

/ d^u^g'^fVauV^n ^ - ^bY = 2^2(1- UaUb) (79) 

(ab) (ab) 

The two terms Eq. and Eq. |79| provide a suitable discretized model 
for a tethered membranes. However, in actual simulations, the even more 
simplified discretization 

F=T.i^a- nf + «^ E (1 - ^i^j) (80) 

is preferred since it is simpler and describes the same universality class (see 



38| for a discussion). Anisotropy may be introduced in this model by ascrib- 



ing distinct bending rigidities to bending across links in different intrinsic 



directions [92 



Self-avoidance can be introduced in this model by imposing that the tri- 
angles that define the discretized surface cannot self- intersect. There are 



other possible discretizations of self-avoidance that we discus in sect. |472 . 

In order to numerically simulate the model Eq. ^ different algorithms 
have been used. A detailed comparison of the performance of each algorithm 
may be found in [|123| . 
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D The crumpling Transition 



The Free energy is given by Eq. [1 
"1 



F(f) 



(81) 



where the dependence on k, is trivially scaled out. The Feynman rules for 
the model are given in Fig.BO. 



m 
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■2 u 

Figure 30: Feynman rules for the model at the crumpling transition 

We need three renormalized constants, namely Z, Z^ and Z^ in order to 
renormalize the theory. We define the renormalized quantities by 

f = (82) 

= M-'Z'^Z~\ , vr = M-'Z'^Z-h . 

Then Eq. becomes 



Fir) 



-{dlvRf + M'Z^u{d^rndprn - ^(5„r^ 



\2\2 



D 



+ 



+ M'z,v{d^fRd''fR) 



(83) 



In order to compute the renormalized couplings, one must compute all 



relevant diagrams at one loop. Those are depicted in fig. |31|. Within dimen- 
sional regularization, diagrams (la) and (lb) are zero, which in turns imply 
that the renormalized constant is Z = 1 at leading order in e, similarly as in 
linear a models. 

Using the integrals in dimensional regularization (see Sect. 0) Diagram 
(2a) gives the result 



_^n*2^jij2 l—{pi- P3P2 -Va + Vi- Pm ■ Ps) + {v^ - ■^)Pi ■ P2P3 ■ Pa 



(84) 
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(1a) 
(2a) 
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Figure 31: Diagrams to consider at one loop 

And diagram (2b) and (2c) may be computed at once, since the result of 
(2c) is just (2b) after interchanging pa p^, so the total result (2a) + (2b) is 



1 1 • • • . ( 61 7 V 

___^^l^2 ^nn ) (^_^2 ' ^3^2 ' P4 + Pi ' Pm ' Pz) 

+ ^"6 +12 + 96^^^ -^^^^-^M 



(85) 



And the result for (2d) and (2e) is just identical, so the total result 
(2d)+(2e) is 



___^n»2^jU2 ) ^ -uv){p^ ■ P2P2 -Pa+Pi- pm ■ Ps) 



13 



u 



+{v + —UV - —)pi ■ P2P3 ■ Pi 



6 



48^ 



(86) 



Adding up all these contributions taking into account the different combi- 
natorial factors (4 the first contribution, 8 the last two ones) and recalling 
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Z = 1, we get 



u 



1 1 Ad 65, 2 4 .\ 

1 1 /21 21 
t^i? + -— ( tt^ur + -^urVr + {Ad + 5)t;^ 



e87r2 Vl6""" ' 2 

The resultant /3-functions are then readily obtained by applying Eq.(^). 



(87) 



E The Flat Phase 

The free energy is given in Eq. ^,and it is given by 
F(u, h)= j c/^x 



The Feynman rules are shown in fig. it is apparent that the in-plane 
phonons couple different from the out-of-plane, which play the role of Gold- 
stone bosons. 
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Figure 32: Feynman rules for the model in the flat phase 
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We apply standard field theory techniques to obtain the RG-quantities. 
Using the Ward identities, the theory can be renormalized using only three 
renormalization constants Z, and Za, corresponding to the wave function 
and the two coupling renormalizations. Renormalized quantities read 



fin = M-'Z'Z-'fi 



u = Z 



(89) 



M-'Z^Z^^X 



Then Eq. ^ becomes 



F(u, h) = / rf^x Z{d^df,hn)^ + 2M'Z^fiRUR^puf + M'Z^K 



a \2 
Ra) 



U 



(90) 

We now compute the renormalized quantities from the leading divergences 
appearing in the Feynman diagrams. The diagrams to consider are given in 
Fig.^ These can be computed using the integrals given in Sect.^. 



(1a) 




(2a) 
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Figure 33: Diagrams to consider at one loop 
The result of diagram (la) is given by 



1 dr 



e 67r2 



/x^(5«/3 - ^) + 3(/.2 + + 2\'f-^^ 



p. 



p. 



P 



(91) 



Diagrams (2b) and (2c) are identically zero, so (2a) is the only additional 
diagram to be computed. The result is 



1 6'^ fijfi + A) 
e87T^ 2/x + A 



10 (p 



2\2 



(92) 
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from Eq. 05] and the definitions in Eq. 



(93) 



Using the previous result in the diagrams (la) whose result is in Eq. ^ we 
obtain 



1 

1 + --r^ifJ^R + G^IrXr + 6\r)/\r 



(94) 



e 247r2 

and we deduce the renormalized couplings 



/i 
A 



1/10 (/i/j + Afi) , 4 \ 2' 

1/10 (/ifi + Ar) c/c , 2 , « X , fix2N 

+ - [-^2^j;^^f^nXR + ^(/i^ + 6/.«A« + 6A^,) 



(95) 



from which the /3-functions trivially follow with the aid of Eq . ([73|) 



F The Self-avoiding phase 

The model has been introduced in Eq. ^ and is given by 

F(r) = ^ / d^x(a„f(x))2 + y rf^xrfV'(r(x) - f(y)) , (96) 

We follow the usual strategy of defining the renormalized quantities by 

f = Z^I'^tr (97) 
h = M'Z.Z^/^R, 

and the renormalized Free energy by 

F(r) = ll rf^x(5,fH(x))2 + M^Z,^ J d^'^dVirRi^) - ^^(y)) • (98) 

The (5-function being non-local adds some technical difficulties to the calcu- 
lation of the renormalized constants Z and Zh- There are different approaches 
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available but we will follow the MOPE (Multilocal-Operator-Product-Expansion), 
which we will just explain in a very simplified version. A rigorous description 
of the method may be found in the literature. 

The idea is to expand the 5-function term in Eq. |98| 



(99) 

with this trick, the delta-function term may be treated as expectation values 
of a Gaussian free theory. This observation alone allows to isolate the poles 
in e. We write the identity 

gifc(r(xi)-r(x2)) _. gjfe(r(xi)-r(x2)) . ^k^G{xi-X2) (100) 

where G{x) is the two point correlator 

- G(x) = -(f(x)f(O)) = j^^^ , (101) 

with Sd being the volume of the D-dimensional sphere. The symbol :: stands 
for normal ordering. A normal ordered operator is non-singular at short- 
distances, so it may be Taylor-expanded 

gifc(r(xi)-r(x2)) _ (l + i(xi -X2)'"(fc9„r) (102) 

-i(xi - X2)"(xi - ^^f{kd^f){kdpf) + ■ ■ .)e^'«(^-^^) . 

To isolate the poles in e we do not need to consider higher order terms as it 
will become clear. If we now integrate over k, we get 

.^(f(x,)-f(x.)) ^ (,,).;.(_^(,^_,^)).,. l (103) 

-\ ( J)"^)^!-^! x,^)l^^^ g/^r(x)g.f(x) + . . . 

= C\^i - X2)l + C"'^(Xi - X2)9^f(x)9„f(x) + ■ ■ • 

where we omit higher dimensional operators in r, which are irrelevant by 
power counting, so since the theory is renormalizable they cannot have simple 
poles in e. Additionally, we have defined x = ^^^^^ . One recognizes in 
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Eq. |1U3| Wilson's Operator product expansion, applied to the non-local delta- 



function operator. 

Following the same technique of splitting the operator into a normal or- 
dered part and a singular part at short-distances, it just takes a little more 
effort to derive the OPE for the product of two delta functions, the result is 

5'^(f(xi) - f(yi))5'^(f(x2) - f(y2)) = C(xi - X2, yi - y2)5'(f(x) - f(y)) + • • • 

(104) 

with ^(xi - X2,yi - y2) = ^^^Y,2^_Gi^^^.l,).G{y^-y2)r/^ - The terms omitted 
are again higher dimensional by power counting so they. The OPE Eq. |103 



and Eq. |104| is all we need to compute the renormalization constants at 



lowest order in e, but the calculation may be pursued to higher orders in e. 
In order to do that, one must identify where poles in e arise. In the previous 
example poles in e appear whenever the internal coordinates (xi and X2 in 
Eq. |103| , xi,X2,y3 and y2 in Eq. |104| ) are pairwise made to coincide. This is 
diagrammatically shown in fig. It is possible to show, that higher poles 
appear in the same way, if more 5-product terms are considered. 




2 



'9 



Figure 34: Diagrammatic expansion to isolate the poles in e within the 
MOPE formalism at lowest non-trivial order. Solid lines represent (5-function 
terms and dashed lines indicate that points inside are taken arbitrarily close. 
Higher orders contributions arise in the same way. 

Let us consider the first delta-function term corresponding to n = 1 in 
the sum Eq. Using Eq. |103| we have 



J d^xdV(ri?(x) - r«(y)) (105) 
Zk J d^xd^y(Ci(x - y) + C"''(x - y)dprR{^)dA{^) + ■■■ 



brM' 
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D 



The first term just provides a renormalization of tlie identity operator, wliicfi 
we can neglect. From 

y|y|>i/Af y D AD s ^ ' ^ ' \V{D/2)) 

(106) 

and we can absorb tlie pole by Z if we define 

From the short distance behavior in the sum Eq. |9^ corresponding to 
77, = 2 we get 

y rf^xf/V'(rR(x) - f^(y)) y t;^zrf^wC(^,«;) , (108) 
where, in order to isolate the pole we can perform the following tricks 

J d^zd^wC{z,w) (109) 



(47r)-'^/24/2(2 - Df' -^f-— ] dz dw 



z w 



V{D/2)) Jo Jo + ^2-D)d/2 

n /a \ 2 D D 



r>/o \ 2 O D 



^ ^ ^ ^ ^ \r(D/2); {2- Df Jx-^+y2<i ^{x + y)-'^l^ 



2 -p/ D ^2 



2-Z)> 



since changing the boundary of integration from a square to a circle does not 
affect the residue of the pole. We finally have 



Z 



h 



1 + £^i(2 - D)-Wi^^_^ 4^ (4vr)-'/' , (110) 

e 2^ ' r(^) \v{D/2)) y ) ^ y ) 
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and the /3-function follows from the definitions Eq. (P7| ) together with Eq. ( |107| ) 
and Eq.(|7|). 



G The mean field solution of the anisotropic 
case 

The free energy has been introduced in Eq. Let us first show the con- 
straints on the couplings so that the Free energy is bounded from below. 

• Uyy > 0: This follows trivially. 

• u'j_j_ = v±± + > : Define A]^ = dav'' (x) then from Eq. ^ we get 

^Tr{AA^r + ^-fiTvAA^r > (-^ + v^^)/2iTrAA^r 
(TrAA^)^ , (111) 



2 

which implies u'^^ > 0. 

v±y > —{u'j_j_UyyY^'^ ! defiuiug b = dyr{x), It is derived from 

{Tr{A^A)f + ^{b^f + v^ybnTrAA^>0. (112) 



Introducing the variables 
and w = ((-D — l)Ci) Cy); the mean field effective potential may be written as 



V{W) = h^-'Ly 



w ■ b + -w ■ A - w 



;ii4) 



In this form, it is easy to find the four minima of Eq. |114| , those are 
1. Crumpled phase: 



Ci = o 



2 _ n Vmrn = (115) 
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2. Flat phase: 



A(D-l) 



tD-It 

4A 



3. ±- Tubule: 



4. Tubule: 



11 

"xx V ■ 

y m.ir. 



u 



±± 



Ci = o 

2 ^ 

3/ Uy 



jD-lj f2 
4 Mj,?/ 



(116) 
(117) 

(118) 



The regions in which each of the four minima prevails depend on the sign 
of A. 

• A > 0: Let us see under which conditions the flat phase may occur. 
We must satisfy the equations 



Uj_J_ty 



< 
< 



'^±y't±y 



(119) 



If v±y > this inequalities can only be satisfied if both t± and tyO have 
the same sign. If they are positive, Eq. |119| imply At± < or Aty < 0, 
which by the assumption A > cannot be satisfied. The fiat phase 
exists for < and < and satisfying Eq. |119| . If ty > and 
t_i_ > then the fiat phase or the tubular cannot exist (see Eq. |117| 
and Eq. |118|) so those are the conditions for the crumpled phase. Any 
other case is a tubular phase, either ±-tubule or ?/-tubule, depending 
on which of the inequalities Eq. |119| is not satisfied. If v±y < it easily 
checked from Eq. |116| that the fiat phase exists as well and the same 
analysis apply. 

A < 0: From inequality Eq. |112| we have v±y > 0. The inequalities are 
now 



U > 



u 



yy 



t± 

v± 

V±y 



U 



(120) 
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Now, in order to have a solution for both inequahties we must have 
^ > which requires A > 0. This proves that the flat phase 
cannot exist. There is then a crumpled phase for ty > and t± > and 
tubular phase when either one of this two conditions are not satisfied. 
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